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Abstract 

We present an analytical calculation of the spectra of CMB anisotropics 
and polarizations generated by relic gravitational waves (RGWs). 
As a substantial extension to the previous studies, three new in- 
gredients are included in this work. Firstly, the analytic Cf^ and 
Cf^ are given; especially the latter can be useful to extract signal 
of RGWs from the observed data in the zero multipole method. 
Secondly, a fitting formula of the decaying factor on small scales is 
given, coming from the visibility function around the photon decou- 
pling. Thirdly, the impacts by the neutrino free-streaming (NFS) is 
examined, a process that occurred in the early universe and leaves 
observable imprints on CMB via RGWs. 

It is found that the analytic Cf^ and Cf^ have profiles agreeing 
with the numeric ones, except that Cf^ in a range / < 10 and the 1^^ 
trough of Cf^ around / ~ 75 have some deviations. With the new 
damping factor, the analytic Cf^^ and Cj^^ match with the numeric 
ones with the maximum errors only ~ 3% up to the first three 
peaks for / < 600, improving the previous studies substantially. 
The correspondence of the positions of peaks of Cf^^ and those 
of RGWs are also demonstrated explicitly. We also find that NFS 
reduces the amplitudes of C^^ by (20% - 35%) for / ~ (100 - 600) 
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and shifts slightly their peaks to smaller angles. Detailed analyses 
show that the zero multipoles Iq, where Ci^^ crosses 0, are shifted 
to larger values by NFS. This shifting effect is as important as those 
caused by different inflation models and different baryon fractions. 

PACS numbers: 98.70.Vc, 04.30.Nk, 98.80.-k, 95.85. Ry 
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1. Introduction 

The observations on CMB [1, 2, 3, 4, 5, 6, 7, 8] are in good agreement with a spatially 
flat universe with nearly scale-invariant spectrum of primordial adiabatic perturbations predicted 
by the inflation model. Generally, two kinds of perturbations of the spacetime metric are of 
interest: density perturbations, i.e. scalar type [9, 10] and relic gravitational wave (RGW), i.e. 
tensorial type [11, 12, 13, 14, 15, 16, 17, 18], respectively. Both perturbations will influence the 
CMB anisotropies and polarizations through the Boltzmann equation for photons. Although the 
contribution from density perturbation is dominant, RGWs may have important contributions [7]. 
In particular, as a special feature, RGWs can give rise to magnetic type of CMB polarizations, 
and this could provide a distinguished way to directly detect RGWs of very long wavelength 
comparable to the Hubble radius ~ I/Hq. In comparison, the usual laser interferometers, such 
as LIGO, probe the intermediate frequency range v = 50 — lOOOHz [19], the waveguide detectors 
probe the high frequency range u — 10^ — lO^Hz [20, 21], and the Gaussian laser beam detectors 
probe in very high frequencies u ~ 10^°Hz [22, 23]. 

The spectra of CMB anisotropies and polarizations generated by RGWs have long been 
computed [24, 25, 26, 27, 28, 29]. In particular, by approximate treatments of the photon 
decoupling, Refs.[30, 31] have derived the analytic expressions of the polarization spectra, Cf^ 
and C^^ , which show explicitly the influences of RGWs, the inflation, the decoupling process, 
the baryons, and the dark energy, etc. But, compared with the numeric computations, these 
analytic Cf^ and C^^ have large errors and are valid only in a limited region of / < 300. This is 
largely due to the damping factor D{k) coming from the visibility function during the decoupling 
is not accurate enough on small scales. Thus, improvements of accuracy and extensions to a 
broader range are certainly desired. We will present a fitting formula of D{k), which substantially 
improves both the accuracy and the region of validity over the previous studies. 

Moreover, in the previous analytic calculations [30, 31], the cross spectrum Cf^ was not given, 
neither was the temperature spectrum C^- Theoretically, the magnetic type of polarization ^ 
can only be generated by RGWs. But the current observed data of C^^ is not yet sufficient to 



confirm the existence of RGWs [4, 5, 7, 6, 8]. On tlie other hand, the cross spectrum Cf^ is 
about two orders higher than Cj^^ , and also contains the contribution from RGWs beside the 
density perturbations. More importantly, WMAP5 has detected Cf^ [7]. This gives rise to the 
possibility to extract RGWs, since the contributions from the scalar and tensorial perturbations 
behave differently. In particular, RGWs can change the value of the multipole /q, where Cf^ 
first crosses 0. Thereby, one can, in principle, tell whether there is contribution from RGWs, and 
give constraints on the scalar/tentor ratio and other cosmological parameters. This so-called 
zero multipole method relies on detailed analysis of Cf^ [32, 33]. Since Cf^ depends on several 
cosmological parameters, analytic results are always helpful in exhibiting their properties. In this 
paper we give the analytic Cf^, as well as Cj"^, due to RGWs. 

As a source to CMB, the RGWs depends on the inflation and on the dark energy [34]. 
Besides, it also depends on physical processes in the radiation-dominated Universe, such as NFS 
[35, 36, 37, 38], the QCD transition, and the annihilation [37, 39, 40]. While the latter two 
processes are effective only on small scales v > 10~^^ Hz and do not appear in the currently 
observed CMB spectra, the former process is effective on large scales with a frequency region 
v ~ (10-^'^ ~ 10-^°) Hz, reducing the amplitude of RGWs by ~ 20% [38]. This in turn will have 
observable effects on the 2"*^ and S^'^ peaks of C^-^ . Thus, we will employ the RGWs spectrum 
modified by NFS [38] to calculate Cf^^, improving the previous calculation in Ref.[31] that did 
not consider the effect of NFS. 

In Section 2 we review briefly the result of RGWs spectrum h{u,r]) with modifications due 
to NFS. In Section 3 we will use this h(i/,r]) to compute the spectra Cf^, Cf^, Cf^ , and 
Cf^. The Basko-Polnarev's method will be used [16, 17, 41, 42]. In the process of the time 
integration, a fitting formula of damping factor D{k) on small scales will introduced, which gives 
a better representation of the visibility function V{ri) during the photon decoupling. Section 4 
examines the influences on C/^^ due to NFS, the spectrum index of inflation, and the fraction 
of baryon, especially, the corresponding shifting of Iq of is investigated. A summary is given 
in Section 5. We use the unit in which c — U — ks — ^ \n this paper. 

2. RGWs Modified by NFS 

The expansion of a spatially flat {Vt^^ + Vt^ + 0.^ = 1) Universe can be described by the 
Robertson-Walker metric 

ds^ = 0^(77) [-dr]^ + (5ij + hij)dx'dx^] , (1) 

where 77 is the conformal time and the small perturbation hij is RGWs and is taken to be traceless 
and transverse (TT gauge) 

h\^0, h%^Q. (2) 



The wave equation of RGWs is 

d^i^d'-hij) = 0. (3) 

By the Fourier decomposition 

h^Av,^)=T.J^,<Jh^\rJ)e^^- (4) 
for each mode k and each polarization a — (+, x), Eq.(3) can be put into the form 

hk + 2-hk + k'^hk = 0, (5) 
a 

where hk — dhk/drj, the polarization index a has been skipped for simplicity. Eq.(5) holds for 
most of the stages from the inflationary to the current accelerating expansion. The explicitly 
form of the scale factor 0(77) are given 

0(77) = ae{r] -rje), Vs < V < V2 (6) 
for the radiation-dominant stage, 

a{v) = ^miv - llmf, V2<V <Ve (7) 
for the matter-dominant stage, and 

a{rj) ^ IhIv - Val"^ , ?7b < ?7 < ?7o (8) 

for the accelerating stage up to the present time r/o, where 7 ~ 1.044 for a dark energy — 0.75, 
and Ih = 1 /Hq with Hq being the Hubble constant. The normalization of airf) is such that 
1^0 — Va\ = 1. where 770 ~ 3-11. The notations in Eqs.(7) and (8) are adopted from that in 
Ref.[34, 38, 39]. In our convention, the conformal time 77 is dimensionless, the scale factor a{vi) 
has the dimension of length. The analytic solution of Eq.(5) and the spectrum are obtained for 
the expanding universe with the consecutive stages: inflationary, reheating, radiation-dominant, 
matter-dominant, and accelerating, respectively in Ref.[34]. 

As also evidenced by the five-year WMAP [5, 7, 6, 8], there exists a cosmic neutrino back- 
ground with the three light species. By the standard scenario of Big-Bang, during the radiation 
stage, from the temperature T ~ 2 MeV up to the beginning of the matter domination, the neu- 
trinos are decoupled from other components and start to freely stream in space. This neutrino 
free-streaming gives rise to an anisotropic portion of the stress tensor, which serves a source 
for RGWs. Consequently, during this period rii,d < f] < "^2, Eq.(5) is modified to the following 
differential-integral equation [35, 36, 37, 38] 

h{ri) + 2-/1,(77) + k^hk{r)) = -24/,(-)2 T h{v')K{kirj - rj'))drj', (9) 



where the kernel of the integral is 



, sinx 3cosx 3smx , , 

Ki^) = -—, 1^ + -::^' (10) 
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and fi, — Pi,/ po is the fractional energy density of neutrinos, whose initial value is f,^{ri = 0) = 
0.40523 for the effective number of species of neutrinos — 3. The term on the right hand 
side of Eq.(9) represents the anisotropic stress tensor due to NFS. 

Eq.(9) has been solved by perturbations, yielding the full analytic solution hk{ri), from the 
inflation up to the present accelerating stage [38], and it has been found that NFS causes a 
damping of hk by ~ 20% in the frequency range 

u ~ (10-^^ 10"^°) Hz. (11) 



Since u is related to the conformal wavenumber k as u = -^k with Hq = 3.24 x 10 ^^h Hz 
being the Hubble frequency, Eq.(ll) corresponds to a range of the conformal wavenumber 



~ (3 X 10\3 X 10^) (12) 

for a Hubble parameter h ~ 0.7. NFS also slightly drags the RGWs spectrum to small scales. This 
dragging effect can be understood by a qualitatively analysis. K{x) in Eq.(lO) has a peak around 
X ~ and K{0) ~ 0.07, and its derivative can be roughly approximated as K'{x) ~ —6{x). 
The integration on the right hand side of Eq. (9) can be integrated by parts 

r hk{v')K{r) - r)')dr)' ~ hk{v){K{0) - 1). (13) 



Then Eq. (9) is approximately reduced to 



a ■ 

hk + 2-hk + 
a 



e - 24/,(i - Km 



hk = 0. (14) 



By comparing Eq.(14) with Eq.(5), one sees that NFS modifies the squared wave number k"^ to 
an effective one ^ 

P = A;' - 24/,(l - i^(0)) <k\ (15) 

If the mode hkirj) without NFS has a peak at A; = kp, then the corresponding mode hk{ri) with 
NFS will have a peak at k{k) = kp, which yields A; ~ [1 + 12/^(1 - X(0)) (ff]kp. The larger 
the kp is, the greater the shifting amount is. This analysis qualitatively explains why NFS slightly 
drags the peaks of RWGs to large k. It is expected that NFS will cause a slight shift of C/^^ to 
large I via RGWs consequently. 

Since the damping range of RGWs is (10~^^, 10~^°) Hz, its lower frequency part just falls into 
the observable domain of C^-^ . Therefore, in calculation of CMB spectra, the RGWs damped 
by NFS should be used as the source. As will be seen in the next section, the mode functions 



hk{r]d) and hk{r]ci) at the photon decoupling time r]d, i.e., z ~ 1100, will appear in the integral 
expressions of the spectra of CMB anisotropies and polarizations. They are plotted in Fig.l. The 
modifications on /ifc(?7d) and hk{r]d) by NFS leave observable imprints in the spectra of CMB. 

As the initial condition, the spectrum of RGWs at the time r]i of the horizon-crossing during 
the inflation is chosen to be [34, 38] 

h{u,r),) = \h,{rj,)\ = A{—f^^-f, (16) 

where kn — Stt is the comoving wave number and corresponds to a physical wave Xh — 
27ra{r]o)/kH — Ih, the constant A is to be fixed by the observed CMB anisotropies, and the 
spectrum index is a parameter determined by inflationary models. The special case of 
An/ — ~2 is the de Sitter expansion of inflation. If the inflationary expansion is driven by a 
scalar field, then the index is related to the so-called slow-roll parameters, 77 and e [43], as 
An/ = — 2 + (?7 — 3e). For demonstration purpose in our context, we allow the parameter An/ 
to take the values > —2. In literature, the RGWs spectrum is often written in terms of A|(/c), 
related to Eq.(16) by h^{v^rii) = 8A^(/c). Without the running index, it is usually assumed to 
have the form [1] [2] [44] 

Al{k)=ATdr. (17) 

Here the tensorial spectrum index ut — '^Pinf + 4, ko is some comoving pivot wavenumber, 
whose corresponding physical wavenumber is ko/a{riH) — 0.002 Mpc~^ [4, 44], and the amplitude 
At = 2.95xlO~^rA{ko) with A{ko) ~ 0.8 as determined by the WMAP observations accordingly 
[2], r being the tensor/scalar ratio. In general, r is model-dependent, and frequency-dependent 
[31, 32]. The value of r has long been an important issue [45, 46, 47, 48, 49]. In our treatment, 
for simplicity, r is only taken as a constant parameter for normalization of RGWs. Currently, 
only observational constraints on r have been given. The 1-Year WMAP gives r < 0.71 [1]. 
The 3- Year WMAP constraint based on the CMB polarization gives r < 2.2 (95% CL) evaluated 
at ko [4], and the full WMAP constraint is r < 0.55 (95% CL) [3, 4]. Recently, the 5-year 
WMAP data improves the upper limit to r < 0.43 (95% CL) [5], and combined with BAO and 
SN gives r < 0.2 (95% CL) [7] [6]. The combination from such observations, as of the Lyman-a 
forest power spectrum from SDSS, 3-year WMAP, supernovae SN, and galaxy clustering, gives 
an upper limit r < 0.22 (95% CL) and r < 0.37 (99.9% CL) [50]. For concreteness, we take 
r ~ 0.37 in our calculation. 

3. Analytical Spectra Of CMB 

In the Basko-Polnarev's method [16, 17], the Boltzmann equation of the CMB photon gas 
for the k-mode is written as a set of two coupled differential equations 

4 + [ikfi + q\^k = hk: (18) 



i3k + [ikii + q\i3k = qGk. (19) 

where Pk represents the linear polarization, = — Pk represents the anisotropy of radiation 
intensity, // = cos^, q is the differential optical depth, and Gk{r]) — ^ /^^^//'[(l + jjI'^YPk — 
1(1 — //'^)^^fe] [31]. Note that the gravitational waves hk in Eq.(18) is the Sachs-Wolfe term 
[51] and plays the role of source to the temperature anisotropies. In the following, we omit the 
subscript k for simplicity of notation. The formal solutions of Eqs.(18) and (19) can be written 
as: 

i(r],ti)= r h{ri')e-<'^^'^'^e'^i'^'^'-'^Ur]', (20) 
Jo 

P(r],fj,)^ r G(r]')q(r]')e-^^''''^'^e''"'^'''-''^drj', (21) 
Jo 

where K,{ri',ri) = jl^ qdrj = k{vi) — Kij]') with Kijf) = K{rio,ri) being the optical depth, such that 
q{rj) = —dK{r]Q,r])/d-r]. To get rid of the angle dependence, { and {3 are usually decomposed in 
terms of the Legendre components 

Uv) = lj'/f^aV,f^)PM, (22) 

f3liv) = lJ'^d^^^3iv,^i)Pli^i), (23) 
where Pi is the Legendre functions. Using the expansion formula 

1=0 

and the ortho-normal relation for Legendre functions, one obtains 

^liVo) = r e-^^^^h{r])Mk{r] - r^,))di^, (24) 
Jo 

PiiVo) = / G{r,)V{r,)3i{k{r, - r,o))dr,, (25) 
Jo 

both being evaluated at the present time t^o, where 

V{n) = q{r))e-<''^ (26) 

is the visibility function for the decoupling. As is known, V{r]) is a narrow function peaked around 
the decoupling time rjd with a width Arja- It phenomenologically describes the details of the 
decoupling process [52, 53, 26]. Ref.[30] uses a single gaussian function to fit V approximately. 
In Ref. [31], as an improvement, the following two pieces of half gaussian function are used 

T/(^,)exp(-(S^), {r,<r„), 
V{V) = { } (27) 



where the decoupling time rjd — 0.0707 corresponding to the redshift — 1100, Arjdi = 0.00639, 
A?7d2 = 0.0117, and {Ar]di + A?7rf2)/2 = Ar]d is the thickness of the decoupling. In absence of 
reionization, the coefficient V{rid) in Eq.(27) will be determined by the normalization 

fVO 

/ Virj)dr]^l. (28) 
Jo 

We have checked that the error between Eq.(27) to the numerically fitted formula given in 
Refs.[52, 26] is very small, ~ 3.9% in the interval r] > rjd. Compared with the single gaussian 
function in Ref.[30], Eq.(27) improves the description of the visibility function by ~ 10% in 
accuracy. Substituting Eq.(27) into Eq.(25), after some treatment of the integration over the 
variable 77 [31], the approximate analytic solution of Pi without reionization has been arrived up 
to the second order of a small l/q^ in the tight coupling limit, 

1 20 

AM = -ln-t'A7jdD{k)h{r]d)ji{k{m-Vo)). (29) 
where h{rid) is the time derivative of RGWs at the decoupling, and 

D{k)^- [e-^''^'^'^^^' + e-^('=^'?d2)^] (30) 

is the Fourier transformation of V{t]) in Eqs.(27) with the parameter c taking values in [0,2]. 
Formally, the occurrence of the damping factor D{k) is due to the i] integration of Eq.(25) of the 

2 

form J^^e~'^ e^^'^dr] = e~~ e"'' dr] since the integrand factor h{r])ji{k{r] — 770)) contains 
a mixture of e*^'' and e^*^''. From view point of physics, D{k) is generically expected [54, 26], 
because the photons diffuse through the baryons around the decoupling and the fluctuations are 
severely damped within the thickness of the surface of the last scattering. Therefore, D{k) is 
very sensitive the thickness Arid- However, as an approximation, the fitting formula Eq.(30) is 
not accurate enough on the small scales and will cause an over-damping of amplitudes of Cf^^ 
for larger This is because, in the afore-mentioned derivation of Eq.(30), other time-dependent 
factors in hk{ri) have been taken as constants during the decoupling. Besides, other processes 
important on small scales were not taken into account [55, 26]. To improve Eq.(30), we adopt 
the following simple fitting formula 

where 6 is a parameter. It will be find that a good fit with c ~ 0.6 and b ~ 0.85, comparing with 
the numerical results. One may even effectively simplify Eq.(31) by the following 

D{k) = e-'^^^''^'^\ (32) 

We find that Eq.(31) and Eq.(32) yield the almost overlapping spectra C^-^ , and the error 
between them is only < 1%. In comparison with the numeric computations, C^^ generated by 
both Eqs.(31) and (32) are much more accurate than those by Eq.(30). 



To evaluate the temperature anisotropies spectrum Cf^, one needs an analytic solution for 
ai. The integrand in (24) contains a factor er'^^^\ which can be treated approximately. Since the 
visibility function V{ri) is a narrow function and can be roughly viewed as a Dirac delta function, 
so by the relation V{r]) — d{e~'^^'^°''^^)/dr], one can treat the factor e~'^^^^ as the step function 



no . 



e-^W ~ / iv<Vd), 

\ 1 {r)d<r)< Tjo). 

Substituting Eq.(33) into Eq.(24) yields 

6(^o) = ^ 
which can be integrated by parts, 

iiiVo) = -i h{Vd)ji{k{r]d - %)) + i / dr]h{r])—ji{k{r] - r]o)), 

JriA ttfl 



(33) 



(34) 



(35) 



where a term containing h{rio) from the upper limit at 770 has been neglected since the amplitude 
h{rio) is about three orders smaller than that of h{rid). The remaining integration term in Eq.(35) 
is small and can be neglected [41], since h{r]) is smaller during the late time {rjd ~ ^^o) and -^ji 
is oscillating functions. Thus one has the following approximate, analytic solution 



6(^0) = -i^h{r]d)ji(k(r]d - m))- 



From Eqs.(29) and (36) follows the temperature anisotropies 

o^iiVo) = -i^ji{k{Vd - Vo)) 



1 20 

h{7]d) + —ln—A7jdD{k)h{7jd) 



(36) 



(37) 



to which both h{r]d) and /i(%) contribute. As our calculation shows, the contribution of h{r]d) 
is about two orders smaller than that of h{r]d). 

In terms of ai and (3i, one calculates C^-^ caused by RGWs [29] straightforwardly. The 
temperature anisotropies 

1 (/ + 2)! ' ' ' ' ^ ^ 2 



k'dk 



ai-2{Vo) 



2a«(^o) 



Stt {l-2)\ 
the electric type of polarization 

1 



+ 



o:i+2ivo) 



{21 - 1) {21 + 1) {21 - 1) {21 + 3) {21 + 1) {21 + 3) 



(3^ 



EE 



167r 



k'dk 



{I + !)(/ + 2)A-2(T/o) ^ 6(/-l)(/ + 2)A(r/o) ^ /(/ - l)A+2(r/o) 



{21 - 1)(2/ + 1) 



(2/-l)(2/ + 3) (2/ + l)(2/ + 3) 



,(39) 



where the second term in the integrand has the coefficient 6{l — 1){1 + 2), different from that in 
Ref.[29], the magnetic type of polarization 



k^dk 



2(/ + 2)A-i(r;o) , 2(/-l)A+i(^o) 



{21 



+ 



{21 



(40) 



and the temperature-polarization cross correlation spectrum 



1 (/ + 2)! 



^ Stt (/-2)! V 167r 



n/-2('/o) 



2n/(//()) 



f>/+2(//o) 



(2i-l)(2i + l) (2/-l)(2Z + 3) (2Z + l)(2/ + 3) 



X 



(/ + 1)(/ + 2)A_2(%) , 6(/-l)(/ + 2)A(7yo) , /(/ - l)A+2(/?o) 



(41) 



(2Z-1)(2Z + 1) (2Z-l)(2i + 3) (2Z + l)(2i + 3) 

where the second term in the integrand has the coefficient 6(/ — 1)(^ + 2), different from that 
in Ref.[29]. Substituting the explicit expressions ai and j3i of Eqs.(29) and (36) into the above 
spectra, one finally has 



^^^^ 1 (/ + 2)! 



Stt {l-2)\ 



1 20 • ^ 

%,) + -ln-Ar^,Z}(A;)M%) , (42) 



EE 



1 



1 



20 



167r V 17 f ^ / k'dkPUKvd - Vo))Ar)lD\k) h(rja) 



(43) 



'^^'^ ^ (]^^"f)7 -^o))Ar^^i52(A;) |%.) 



(44) 



_ (1+2) 

136v^7r 3 ^ (/ - 2) 



1 

2 

+ 



1 20 

-Kvd) -—In— ArjdD{k)h{ria) 
1 20 

-/.*(%)-- In -A%D(A;)/i*(%) 



^ 1 edkPTi(k(rid - Vo))PEi(k(rid - Vo)) 
h{r]AAr]dD{k), 



where the projection factors are defined as [42]: 

ji-2{x) , 2ji{x) 



Pti{x) = 



+ 



+ 



(2/-l)(2/+l) (2/-l)(2/ + 3) (2/ + l)(2/ + 3) x 



(45) 



(46) 



(2Z-1)(2Z + 1)- 



{21 - i){2i + ?>y 



(2Z + l)(2Z + 3)' 



X 



, , 2(/ + 2) , , 2(/ - 1) , , , , / - 1 . ^ 



(47) 



(48) 



With h{rid) and /i(77d) given from the last section and D{k) from Eq.(31), we compute Cj^-^ 
and plot them in Fig. 2, where the following values of respective parameters are taken: the 
inflationary index Pinf — —2.02, the dark energy JIa = 0.75, the baryon density = 0.045, 
the neutrino species Ni, = 3, the tensor/scalar ratio r = 0.37, c = 0.6, and h — 0.85. For 
comparison, in Fig. 2 the numerical results from CAMB [25] are also plotted. It is seen that the 



analytic Cf^ and Cj^^ agree very well with the numerical ones for the range / < 600 covering 
the first three peaks, and the error is only ~ 2%. Comparing with the previous analytic evaluation 
in Refs.[30, 31], our result not only extend the range of validity from / < 300 to / < 600, but 
also improves accuracy substantially. 

For the spectra Cf^ and Cf^, the profiles of our analytic result also agree with the numerical 
ones fairly well, except that C-f^ in a range / < 10 and the trough of Cf^ around I ~ 75 have 
some deviations. For the purpose of extracting RGWs, more important is Cf^ , whose amplitude 
at the 1** trough has a maximum deviation ~ 20% from that of the numerical CAMB. This is 
due to the approximation of the temperature anisotropies ^;(77o) 'n Eq.(36), which is not accurate 
enough for very large scales. 

The profiles of C/^^ are largely determined by those of h{rid) and h{rid), especially, the peaks 
and troughs of C/^^ correspond to those of RGWs. The integrands for Cf^^ in Eqs. (42), 
(43), (44), and (45) contain the respective projection factors, Pti, Peu Pbu which are made of 
the spherical Bessel's functions. Since ji{x) is rather sharply peaked around x ~ Z for large I, 
consequently, the projection factors as functions of k are peaked around 

K'Ho - Vd) ~ /c77o - ^- (49) 

Therefore, C^^ as integrations over k will receive major contributions from the integration 
domain k ~ l/r]o [31]: 

C^roc|M%)lL/.o' (50) 

Cr,Cr^\h{ria)\l^^^^^D\k), (51) 

oc h{r)d)h{rid)k^i/r,oD{k). (52) 

By Eq.(50), the locations of the peaks of Cf^ is mainly determined by \h{r]d)\'^. Indeed, the 
right panel in Fig. 3 shows that the peaks and troughs of Cf^ correspond to those of \h{rid)\'^. 
Similarly, by Eq.(51), the locations of the peaks of Cf^ and C^^ correspond to those of |/i(%)P> 
shown in the left panel in Fig. 3. The similar correspondence, as revealed by Eq.(52), of the 
locations of the peaks of C^^ to those of h(r]d)h(r]d)k is also confirmed, but the graph is not 
presented though in order to save room. Besides, Cf^, C^^ and Cf^^ also depend on the 
damping factor D(k), leading to strong damping of Cf^, Cj^^ and Cf^ at large /, as shown 
in Fig. 2 and Fig. 3. Overall, our analytic formulation yields a good approximation of C/^^ in 
comparison with the numerical results. 

We use the data from 5-year WMAP [7, 8] to limit the B-mode polarization Cj^^ generated 
by RGWs in Fig. 4. It is seen that the current observational data can only put a rather loose 
limit on RGWs. For r = 0.37 the amplitude of the predicted spectrum is about 2 orders below 
the upper limit by WMAP5. Improvements on the limit, or possible direct detections of C^^ 



are expected from more sensitive polarization measurements by upcoming experiments, such as 
Clover [57], EBEX [58], QUIET [59], Spider [60], and Planck [61]. 

4. Influences by NFS, Inflation, and Baryons 
The NFS 

Let us analyze the effect of NFS on the spectra C^-^ . To demonstrate this, the spectra 
Ci^-^ with and without NFS are plotted in Fig. 5. The I < 100 portion of the spectra are not 
much affected by NFS, only on the scales of / > 200, are the spectra modified effectively. The 
reduction of amplitudes of Cf^, Cf^ and Cf^ by NFS is noticeable only starting from the 
second peak. For instance, the third peak of Cf^ is reduced by ~ 25% and the fourth peak by 
~ 35%. Similar modifications also occur in the spectra C^^, Cf^, and Cf^. These features 
of modifications can be understood as follows. As shown in Eq.(12), the damping of RGWs is 
effective for the conformal wavenumbers k > 30, which, by the relation in Eq.(49), means that 
only those portion with I > 100 of CMB spectra will be affected by NFS. Besides, Figs. 5 also 
shows that NFS causes a slight shift of the peak locations of C^-^ to larger I, a feature to be 
expected, since NFS shifts the peaks of h{rid) and h{'qd) slightly to large k, the peaks of C^-^ 
will be accordingly shifted to larger / by Eqs.(49) - (52). Given the current precision level of 
observations on CMB, these small modifications caused by NFS will be difficult to detect at the 
moment. 

As mentioned in the Introduction, the cross spectrum Cf^ can be useful in revealing the 
presence of RGWs in the zero-multipole method [32, 33, 41, 42, 56]. The 5-year WMAP [7, 8] 
has given the observed Cf^, which is negative (anti-correlation) in a range / ~ (50,220). 
Theoretically, it is a combination of contributions of the density perturbations and the RGWs as 
well. To search for the evidence of RGWs, one needs to disentangle the contribution of RGWs 
from the total. The inclusion of NFS into the calculation will cause a shifting of the position 
of the peaks of Cf^ to larger values of I, and AZ tends to increase with I. For instance. Fig. 5 
shows that, without NFS, Cf^ < for / < 136 and Cf^ > for Z ~ (137, 179). When NFS is 
included, Cf^ < for Z < 136, Cf^ > for Z ~ (137, 183). In this low I region the shifting 
due to NFS is small AZ < 4. But, in the large I region the shifting is large, say, around I ~ 500, it 
is AZ ~ 10. This analysis tells us that the zero multipole around Iq ~ 50 is not strongly affected 
by NFS. However, if we look at the the second zero multipole I ~ 220, at which Cf^ crosses 
once again and turns positive, the influence by NFS is rather strong, AZ ~ 5. More accurate 
observations of Cf^ and detailed analysis are needed, before a definite conclusion can be drawn 
on the existence of RGWs. 

The inflation 

The CMB spectra generated by RGWs depend very sensitively on the initial spectrum h{v, r]i) 
during the inflationary stage. For the power-law form of h{v,r]i) in Eq. (16), Cf-^ depend on 



both the amplitude A and the index Fig. 6 shows Q for the cases of = —1.8 and 
—2.02 with NFS being taken into account. A larger index P^f yields higher amplitudes of Cf^ 
and in the whole range of /, agreeing with the previous result [34], and higher amplitude 
of Cf"^ for the range / > 20. In the zero multipole method, one is more interested in Cf^ in 
the narrow range / ~ (40,60), in which the first zero multipole Iq should appears. Firstly, as 
Fig. 6 shows, Cf^ is negative in this range, and moreover, a larger index (3inf{= —1.8) yields a 
steeper, down slope of Cf^ of negative amplitude. With other parameters being fixed, a larger 
index Pi^f tends to shift the value of the zero multipole Iq of Cf^ to larger I. For instance, our 
calculation shows that, relative to the WMAP-preferred Pinf — —2.02 case, the exact de Sitter 
Pinf — —2 case shifts lo to a larger value by Al ~ 1, and the less-preferred case Pinf — —1.8 
shifts lo by Al ~ 8. 

Notice that, for Cf^^ around I ~ 20, the two curves for P^f — —1.8 and for Pinf — —2.02 
intercept. A similar interception also occurs for Cf^ as well. This behavior can be understood 
as the following. The initial spectrum of RGWs in Eq. (16) contains a factor {^Y'^^^"^ with 
the comoving pivot wavenumber ko ~ 6h~^. On large scales k < ko, one has {-^Y'^^'"^ < 1 for 
p-^j: = -1.8, and (^)^+^*"^ > 1 for pinf = -2.02. By the relation in Eq.(49), the corresponding 
pivot multipole is I ~ krjo ~ 20. Thus, in the region of / < 20, Cf^ and Cf^ have lower 
amplitude for Pinf = —1.8, and higher amplitude for Pinf = —2.02. 

The baryon density 

The wave equation (5) of RGWs is not explicitly coupled with the baryons. As a result, 
hk{r]) and hk{r]) are not very sensitive to the baryons. The impact on CMB by the baryons are 
mainly through the Thompson scattering terms, q^k, QPk, and qGk, in Boltzmann's equation 
of photons (18) and (19). During the evolution of CMB, the photon decoupling process is 
particularly important, which depends sensitively on the baryon component. The fitting formula 
of the visibility function V{r]) given in Ref.[26] contains explicitly the baryon fraction A larger 
ilb yields a larger decoupling time rjd and and a smaller decoupling width At^^ [31]. Moreover, in 
Eqs. (42) through (45), the integrands contain h{rjd), h{rid), and ArjdD{k), which are functions 
of rjd and Arjd- We plot C^-^ for Jl^ = 0.045 and 0.1 in Fig. 7. It is seen that the amplitudes 
of Cf ^ and Cf^ with = 0.045 is slightly higher than those of Qb = 0.1. Thus, a larger 
gives a lower amplitude of Cf^ and C^^, agreeing with the previous calculations [31]. As a new 
result of this paper. Fig. 7 also shows that a smaller 0,^ yields a higher amplitude of Cf^ and 
shifts the value of the zero multipole Iq to large I. For instance, Qf, — 0.045 shifts Iq to a large 
value by AZ ~ 2 relative to the = 0.1 case. Besides, Cf^^ is less sensitive to the value of 0,b 
than the other three spectra. 

5. Summary 

In this paper we have presented the approximate, analytical formulae of the four CMB spectra 



generated by RGWs. This has been motivated by an attempt to extract signals of RGWs possibly 
already contained in C/^"^, especially in the magnetic polarization spectrum Cj^^ and the cross 
spectrum Cf^ . 

In our calculation, a fitting formula of the exponentially damping factor D{k), in Eqs.(31) or 
(32), has been introduced to describe the decoupling process effectively. The resulting analytic 
spectra Cf^ and Cj^^ agree quite well with the numerical ones from CAMB on large scales for 
the first three peaks for I < 600, and the error is only ~ 3%. This improves substantially both the 
precision and the range of validity in comparison with the previous analytic studies. The spectra 
Cf^ and Cf^ are first analytically computed in this paper. They have overall profiles agreeing 
with the numerical ones, but their amplitudes have certain deviations due to the approximation 
adopted in Eq.(36). More relevant to us is Cf^, whose amplitude of the 1** trough at Z ~ 75 
has a maximum deviation ~ 20%. An analytic formulation of ^(770) better than Eq.(36) should 
be aimed at in future work. 

For the Sachs-Wolfe term in the Boltzmann's equation for photons, we have included the 
damping effect of NFS on the RGWs hk{ri) as the source. As is expected, NFS appreciably 
reduces the amplitudes of C/^^ for large I > 100 and, at the same time, shifts slightly the 
locations of the peaks to large I. Thus, in the zero multiple method by examining the positions 
where Cj^ crosses 0, the shifting due to NFS effect should be taken into account for a complete 
analysis. 

We have also demonstrated the influences on Cf^^ by the tensorial spectrum index fimj of 
the inflation and the baryon fraction Vti,. It is found that a larger leads to a higher amplitude 
of CMB spectra, whereas the larger Vti, gives a lower one. Both of them shift the locations of the 
peaks of C/^^. In regards to the shifting of the zero multipoles /q of , NFS is as important 
as the inflation and the baryons and should be included in any comprehensive study. 
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Figure 1: The RGWs hk{rjd) and hk{rjd) at the decouphng. NFS reduces the amphtudes and 
shifts the peaks to larger k. 



Figure 2: The analytic spectra C^-^ generated by RGWs are compared with the numeric 
ones from CAMB [25]. Here the decay factor D{k) in Eq.(31) has been used. 



Figure 3: The left panel: the locations of peaks Cj^^ and C/^^ correspond to that of \hk{rjd)\'^. 
The right panel: the locations of peaks of Cf'^ correspond to that of \hk{rid)\'^. Here \hk{rid)\'^ 
and \hk{rid)\'^ have been plotted with the variable kriQ, which is ~ Z by Eq.(49). 



Figure 4: The predicted Cj^^ is well below the constraint of the 5-year WMAP data [7, 8]. 
Here the tensor/scalar ratio r = 0.37 is taken in computation. 



Figure 5: The NFS modifications on Cf^^ are demonstrated. For I < 600, the amphtudes 
are reduced by up to ~ 35% and the peaks are shifted to larger I by NFS. 



Fierure 6: Cf"" are sensitive to the inflation index /3j„j of RGWs. Two case are plotted for 
An/ = — 1.8 and —2.02. A larger yields higher amplitudes of Cf^ and Cf^. 



Figure 7: The baryon fraction fife affects C^^ . Two cases are given for = 0.1 and 0.045. 
A larger Qf, yields lower amplitudes of Cf^, Cf^, and Cf^. 
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